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Introduction 

The  ability  of  near-infrared  (NIR)  light-based  techniques  to  noninvasively  image 
and  analyze  tissue  structure  and  function  promises  their  great  potential  for  detection 
and  diagnosis  of  breast  cancer.  Optical  diagnostic  techniques  allow  us  to  not  only 
enhance  the  existing  capabilities,  but  to  eliminate  the  need  for  physical  biopsies.  In 
addition,  optical  imaging  is  inexpensive  and  portable,  which  indicates  that  optical 
imaging  could  be  an  ideal  candidate  for  routine  breast  screening.  However,  since  the 
scattering  properties  of  tissues  convolute  re-emitted  NIR  signals,  the  extraction  of 
pertinent  information  continues  to  remain  elusive.  An  understanding  of  light  propagation 
and  light-tissue  interaction  is  required  before  the  optical  technologies  can  substantially 
impact  diagnostic  medicine. 

Research  efforts  in  the  Biomedical  Optics  Laboratory  at  Clemson  University  are 
focused  on  the  biophysics  of  light  propagation  and  light-tissue  interaction  in  order  to 
engineer  appropriate  approaches  for  noninvasive  breast  imaging  and  spectroscopy. 
Specifically,  we  are  developing  indirect  optical/fluorescence  approaches  using  photon 
migration  measurements  in  the  continuous-wave  and  frequency  domains.  These 
indirect  optical/fluorescence  approaches  or  image  reconstructions  are  computationally 
based  on  the  powerful  finite  element  methods.  A  CCD-based  optical  spectroscopic 
imaging  system  is  already  operational  in  our  laboratory  for  continuous-wave 
tomographic  photon  migration  measurements,  while  a  frequency-domain  system  is  still 
under  construction.  Using  these  optical  systems  coupled  with  our  finite  element  based 
reconstruction  algorithms,  we  will  be  able  to  extract  spatial/spectroscopic  maps  of  tissue 
optical  properties,  lifetime  and/or  yield  of  endogenous  and  exogenous  fluorescent 
probes.  Since  metabolic  tissue  states  can  be  identified  by  our  approaches,  diagnostic 
information  is  also  obtained  in  addition  to  detection  of  tumor. 

This  Career  Development  application  for  support  of  Dr.  Huabei  Jiang  will  facilitate 
the  establishment/continuation  of  these  research  activities.  Interdisciplinary  interactions 
with  the  Greenville  Hospital  System  (Greenville,  SC)  will  be  enhanced,  which  insures 
the  direction  of  research  towards  a  clinically  pertinent  and  feasible  system. 

Body 


This  report  describes  work  accomplished  during  the  second  year  of  a  proposed 
four-year  study.  This  Career  Award  supports  Dr.  Jiang’s  research  on  optical  and 
fluorescence  imaging  using  both  continuous-wave  and  frequency-domain 
measurements.  The  focus  of  the  proposed  work  in  Year  2  is  the  continued  evaluation 
and  optimization  of  some  enhancing  schemes  in  our  existing  2D  reconstruction  codes 
using  tissue  phantom  experiments;  continued  phantom  studies  of  image  reconstructions 
with  dye-free  or  dye-laden  phantom  background;  improvement  of  the  imaging  system 
that  was  built  during  the  first  year  of  this  project;  in  vivo  animal  imaging  studies  using 
frequency-domain  measurements. 

Software  work:  We  have  continually  evaluated  and  optimized  the  important  image 
enhancing  schemes  that  were  implemented  during  the  first  year  of  the  project. 
Extensive  tissue  phantom  experiments  have  been  conducted  to  complete  these  studies. 
The  results  from  these  studies  have  been  published  in  the  peer-reviewed  literature  (see 
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a  manuscript  entitled  “Quantitative  optical  image  reconstruction  of  turbid  media  by  using 
direct-current  measurements”  provided  in  the  Appendix  to  the  Summary  Report). 

Hardware  Development:  We  have  attempted  to  improve  the  multi-channel  frequency- 
domain  imaging  system  we  constructed  in  Year  1  in  terms  of  the  signal-to-noise  ratio. 
What  we  have  tried  was  to  use  a  thermo-electric  cooling  system  for  the  photomultiplier 
tubes  (PMTs).  This  has  provided  20-30%  noise  reduction  relative  to  the  original  imaging 
system.  We  plan  to  continually  improve  the  system  in  other  areas  such  as  the  use  of 
low  loss  silica  optic  fibers  (the  current  system  uses  glass  optic  fibers  which  normally 
give  30-40%  more  losses  than  the  silica  fibers). 

Phantom  Experiments:  Using  our  frequency-domain  imaging  system,  we  have 
conducted  extensive  phantom  experiments  for  fluorescence  lifetime  reconstruction  with 
dye-free  background  and  dye-laden  background  which  was  proposed  to  take  place 
during  the  third  year  of  the  project.  Our  experimental  setup  used  was  the  automated 
multi-channel  frequency-domain  system  mentioned  above.  The  system  employed  a 
radio-frequency  intensity-modulated  near-infrared  beam.  The  laser  beam  at  785  nm  was 
sent  to  the  phantom  by  16  fiber  optic  bundles  coupled  with  a  high  precision  moving 
stage.  The  diffused  radiation  was  received  by  another  16  channel  fiber  optic  bundles 
and  delivered  to  a  thermo-electric  cooled  PMT.  A  second  PMT  was  used  to  record  the 
reference  signal.  These  PMTs  were  supplied  at  a  radio-frequency  modulated  current 
with  0.1-1  KHz  shift.  The  intermediary  frequency  signal  obtained  from  the  PMTs  was 
processed  using  a  National  Instruments  board.  To  increase  the  dynamic  range  of 
intensities  that  the  PMT  can  detect,  an  automated  filter  wheel  with  pre-calibrated  neutral 
density  filters  was  added  to  the  system.  For  every  source  position,  16  measurements  for 
each  detector  were  made,  taking  alternatively  100  ms  samples  for  sample  and 
reference  signals.  Fluorescence  signals  were  obtained  through  an  830nm  interference 
filter  placed  in  front  of  the  detection  PMT.  ac  intensity  and  phase  shift  between 
reference  and  sample  signals  were  obtained  using  FFT  Labview  routines.  The  total  data 
collection  time  for  256  measurements  was  8  minutes. 

Solid  phantom  was  used  to  mimic  the  human  tissue.  It  was  made  of  agar, 
Intralipid,  black  ink  and  fluorescent  dyes.  The  absorption  and  the  reduced  scattering 
coefficients  are  linear  with  the  ink  and  Intralipid  concentrations,  respectively. 
Micromolar  ICG  and  DTTCI  dyes  were  added  in  the  tissue  phantom  to  provide  the 
fluorescence  contrast.  Agar  is  used  to  make  the  phantom  solid.  The  solid  phantom 
consisted  of  a  cylindrical  background  and  a  cylindrical  heterogeneity.  The  absorption 
peaks  of  ICG  and  DTTCI  are  764  nm  and  780  nm,  respectively.  And  the  fluorescent 
emission  peaks  of  them  are  803  nm  and  830  nm  respectively.  Their  lifetimes  in  water 
were  measured  to  be  0.56  ns  and  1.18  ns,  respectively.  A  laser  diode  of  wavelength 
785  nm  was  used  to  excite  both  dyes,  and  the  emission  light  of  wavelength  at  830  nm 
were  detected  for  both  of  them  through  an  interference  filter  of  center  wavelength  830 
nm  with  lOnm  bandwidth.  The  image  results  have  been  presented  in  the  second  and 
third  manuscripts  provided  in  the  Appendix  to  this  Summary  Report. 

Initial  In  vivo  Experiments:  We  have  tested  our  overall  reconstruction  approach  using 
in  vivo  animals  which  we  originally  planned  to  pursue  after  the  funded  four-year  study. 
While  we  may  be  able  to  conduct  another  in  vivo  experiment  within  the  course  of  this 
project,  this  initial  animal  study  has  surely  provided  us  an  excellent  opportunity  to  test 
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the  feasibility  of  the  proposed  research  approach  for  breast  imaging.  In  the  in  vivo 
experiments,  Wistar  furth  female  rats  with  MT/W9a-B  (mammary  carcinosarcoma  )  cell 
line  developed  in  their  right  udder  were  used  as  the  animal  models.  Animals  were 
studied  when  the  tumor  sizes  were  around  10-30  mm  diameters.  Animals  were  initially 
anesthetized  with  2%  isoflurane.  The  femoral  vein  was  cannulated  for  the  administration 
of  intravenous  drugs,  including  the  dye  ICG.  Anesthesia  was  maintained  during  the 
experiment.  The  animals  were  placed  in  a  holder  with  their  tumor  positions  placed  into 
the  measuring  ring.  1.5  mg/Kg  bodyweight  of  ICG  dye  was  injected  into  the  animal. 
Imaging  studies  were  carried  out  after  2  minutes  of  the  dye  injection.  Successful  in  vivo 
image  results  have  been  reported  in  the  third  manuscript  provided  in  the  Appendix. 

Key  Research  Accomplishments 

1.  We  have  continually  evaluated  and  optimized  the  image  enhancing  schemes 
developed  in  Year  1  using  tissue  phantom  experiments. 

2.  We  have  improved  the  multi-channel  frequency-domain  imaging  system  that  was 
constructed  in  Year  1 . 

3.  We  have  conducted  extensive  phantom  experiments  for  fluorescence  lifetime 
imaging.  The  successful  experiments  have  confirmed  the  imaging  capability  of  our 
reconstruction  software.  We  have  also  performed  successful  fluorescence  phantom 
studies  that  were  proposed  to  occur  in  Year  3  of  the  project. 

4.  We  have  for  the  first  time  obtained  fluorescence  lifetime  images  from  in  vivo  animal 
measurements  which  were  even  not  proposed  to  attempt  in  the  original  proposal. 

Reportable  Outcomes  (see  the  Appendix  to  this  Summary  Report) 

Conclusions 

We  have  made  a  significant  progress  that  has  exceeded  the  statement  of  work 
proposed  for  Year  2  of  this  project.  Given  the  successful  first  year,  we  will  be  able  to 
fulfil  or  exceed  the  work  statement  for  Year  3,  and  we  will  carry  these  successful  works 
over  into  the  final  year  of  this  project. 
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Quantitative  optical  image  reconstruction  of  turbid 
media  by  use  of  direct-current  measurements 


Nicusor  Iftimia  and  Huabei  Jiang 


We  present  a  detailed  experimental  study  concerning  quantitative  optical  property  reconstruction  of 
heterogeneous  turbid  media  by  use  of  absolute  dc  data  only.  We  performed  experiments  by  usmg 
tissuelike  phantoms  in  both  single-target  and  multitarget  configurations  in  which  variations  m  target 
size  and  optical  contrast  with  the  background  were  explored.  Our  results  show  that  both  scattering  and 
absorption  images  can  be  reconstructed  quantitatively  by  use  of  dc  dafcnmly,  whereas  it  was  impossible 
to  obtain  such  quantitative  information  in  previously  reported  studies.  We  believe  that  this  improve¬ 
ment  is  primarily  a  result. of  the  realization  of  a  novel  data  preprocessing/optimization  scheme  for 
accurately  determining  several  critical  parameters  needed  for  reconstruction.  The  use  of  this  data 
-preprocessing/optimization  scheme  also  eliminates  the  calibration  reference  measurement .  previously 

_  required  for  reconstruction.  Experimental  corrfircnation“of  this  scheme  is  given  in  detail.  ©  2000 

Optical  Society  of  America 

OCIS  codes:  170.3010,  170.3830,  170.3890,  170.6960. 


1.  Introduction 

Near-infrared  diffusive  light  can  be  used  to  image 
thick  tissues  such  as  breast  and  brain,  as  shown  by 
recent  significant  progress  in  the  area  of  optical 
imaging.1-3  Among- various  optical  methods  that 
are  being  developed,  reconstruction-based  optical  im¬ 
aging  has  received  particular  attention  because  it  can 
extract  both  structural  and  functional  information 
from  tissue.  To  date  various  reconstruction  algo¬ 
rithms  have  been  tested  and  evaluated  by  using  con¬ 
siderable  experimental  data  from  different  hardware 
systems  including  cw,  frequency,  and  time 
domain.4-17  Although  significant  progress  has  been 
made,  a  critical  need  to  explore  and  improve 
reconstruction-based  optical  imaging  performance  in 
experimental  conditions  remains. 

We  are  interested  in  frequency-domain  systems  in 
part  because  of  the  presence  of  a  measurable  phase 
component  that  can  be  exploited  during  image  forma¬ 
tion  as  well  as  cw  or  dc  systems  because  of  the  added 
simplicity  of  recording  only  intensity  data.  When 
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powerful  finite  element-based  reconstruction^  algo¬ 
rithms  were  used,18  quantitative  reconstruction  of 
both  absorption  and  scattering  images  has  been 
achieved  from  frequency-domain  measurements!7*19"21 
From  de-data,  simultaneous  recovery  of  both  absorp¬ 
tion  and  scattering  images  has  been  obtained,  but 
^quantitative  reconstruction  in  terms  of  optical  prop¬ 
erty  values  was  impossible.22*23  In  these  early  studies 
we  used  a  calibration  procedure  based  on  a  homoge¬ 
neous  medium  measurement  to  calibrate  the  boundaiy 
conditions  (BC’s)  and  the  source  term.  Although  it 
was  straightforward  to  do  so,  a  trial-and-error  process 
was  needed  to  complete  this  calibration  procedure. 
We  found  that  the  amplitude  of  the  source  largely  de¬ 
termined  the  computed  photon  density  amplitude  and 
that  BC’s  could  affect  the  computed  overall  light  dis¬ 
tribution  considerably,  which  suggests  that  this  time- 
consuming  procedure  may  introduce  significant  errors 
into  the  reconstruction  if  not  carefully  done.  More 
'  important,  it  may  not  be  practical  to  perform  such  a 
calibration  measurement  in  an  in  vivo  setting.  .  We 
also  believe  that  the  inaccuracy  caused  by  this  calibra¬ 
tion  procedure  made  it  impossible  to  obtain  quantita¬ 
tive  reconstruction  from  dc  data,  as  confirmed  by  this 
study. 

Although  hampered  by  the  requirement  of  a  cali¬ 
bration  measurement  in  the  past  few  years,  we  have 
recently  realized  a  novel  data  preprocessing/ optimi¬ 
zation  scheme  that  allows  for  accurately  determining 
the  critical  parameters  needed  for  reconstruction  in¬ 
cluding  the  BC  coefficient,  the  excitation  source  term, 
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and  the  initial  estimates  of  optical  properties  without 
reference  or  calibration  measurements.  In  this  pa¬ 
per  we  present  in  detail  an  experimental  confirma¬ 
tion  of  this  new  scheme  by  using  dc  measurements. 
In  addition,  we  demonstrate  that  the  use  of  this 
scheme  allows  for  quantitative- recovery  of  both  ab¬ 
sorption  and  scattering  images  from  dc-only  mea¬ 
surements.  Experiments  with  tissuelike  phantoms 
in  both  single-target  and  multitarget  configurations 
have  been  performed  where  variations  in  target  size 
and  optical  contrast  with  the  background  have  been 
explored. 

2.  Methods  and  Materials 


A.  Data  Preprocessing/Optimization  Scheme 


We  have  used  in  our  study  a  finite-element-based 
algorithm  for  image  reconstruction  that  has  been  pre¬ 
viously  described  in  detail.18'19  The  algorithm  uses 
a  regularized  Newton  method  to  update-  an  initial 
(guess)  optical  property  distribution  iteratively  to 
minimize  an  object  function  composed  of  a  weighted 
sum  of  the-squared  difference  between  computed  and 
measured  data.  To  illustrate  the  data  preprocess-^ 
ing/ optimization-scheme,  here  we  briefly  outline  our 
reconstruction  algorithm 


Since  in  this  paper  we  focused  on  reconstruction  by 
using  dc  data,  our  image  reconstruction  algorithm  is 
based  on  the  steady-state  diffusion  equation,  which 
can  be  stated  as 


V-D(r)V<E>(r)  -  na(r)<P(r)  =  -S(r),  Q) 

where  <t>(r)  is  the  photon  density,  |j.a(r)  is  the  absorp¬ 
tion  coefficient,  and  D(r)  isjhe  diffusion  coefficient, 
which-can  be  written  as  D{r)  =  l/3[|ia(r)  4-  p.s'(r)],’ 
where  p.s'(r)  is  the  reduced  scattering  coefficient! 
S(r)  is  the  source  term.  In  this  study,  a  point  source, 

S  =  S08(r  -  r0),-is  used  in  which  S0  is  the  source 
strength  and  S(r  —  r 0)  is  the  Dirac  delta  function  for 
a  source  at  r0. 

It  has  been  shown  that  the  following  type  EH  BC’s 
can  provide  the  most  accurate  solution  of  Eq.  (I)19.2*; 

-DV<t  •  n  =  a<J>,  (2) 

where  h  is  the  unit  normal  vector  for  the  boundary 
surface  and  a  is  a  coefficient  related  to  the  internal 

reflection  at  the  boundary.  _ 

Making  use  of  a  finite  element  discretization,  we 
can  obtain  a  discrete  matrix  for  Eq.  (1)  and  realize 
other  derived  matrix  relations  through  differentia¬ 
tion,  which  lead  to  a  set  of  equations  capable  of  an 
inverse  problem  solution7'18'19: 


[Am = {b},  (3) . 


[A] 


(4) 


(Sr3  +  A2) AX  -  3Wm)  -  <1>,CI],  (5) 

^here  the  elements  of  matrix  [A]  are  ay  =  (-DV6-  • 
■  <'  ~  M-<A  <!>/>,  where  ( )  indicates  integration  over  tlie 


problem  domain,  (f>(  and  ct >j  are  locally  spatially  vary¬ 
ing  Lagrangian  basis  functions  at  nodes  i  and  j ,  re¬ 
spectively;  x  expresses  D  or  pa;  3  is  the  -Jacobian 
matrix  that  should  be  formed  by  d$/dx  at  the  bound¬ 
ary  measurement  sites;  Ax  =  (A Dv  AD2,  .  . . ,  A DN, 
Ap-a.n  Ap.a  2, .  . . ,  A is  the  update  vector  for  the 
optical  property  profiles,  where  N is  the  total  number 
of  nodes  in  the  finite  element  mesh  used-  = 

‘  H 


<Vm;]_  and  <E>(C)  =’[$1^  ^,c> 


Vc)],  where  4>/m)  and  T/c\  respectively,  are  mea¬ 
sured  and  calculated  data  for  i  =  1,  2,  .  . . ,  M  bound¬ 
ary  locations.  Note  that  to  estimate  D  and  |xa 
spatially,  we  need  to  expand  these  quantities  in  a 
similar  manner  to  $  as  a  finite  sum  of  unknown 
coefficients  multiplied  by  the  locally  defined  Lagrang¬ 
ian  basis  functions. 


In  optical  image  reconstruction,  the  goal  is  to  up¬ 
date  the  D  and  p.a  distributions  through  the  solution 
of  Eqs.  (3)— (5)  so  that  a  weighted  sum  of  the  squared 
difference  between  computed  and  measured  data  can 
be  minimized^  Note  that  four  critical  parameters 
(the  BC  coefficient  a,  the  source  strength  S0,  and  the 
initialjpiesses  of  D  and  |ra)  must  be  accurately  de¬ 
termined  to  obtain  high-quality  reconstructions!-  a 
is  related  to  reflection  at  the  tissue-air-  detector  in¬ 
terface;  it  is  impossible  to  obtain  an  analytic  expres¬ 
sion  of  a  for  a  realistic  situation,  particularlyfor  an  in 
vivo  sitting.  Although  in  principle  S0  can  be  mea- 
sured^in  a  realistic  situation  it  is  time-consuming 
and  inaccurate  to  do  so.  Although  it  is  usually  pos¬ 
sible  to  make  a  good  initial  estimate  of  D  and  for 
laboratory  experiments,  it  will  be  difficult  to  make 
such  estimates  in  clinical  situations.  As  indicated  in 
Section  l_our  early  research  in  this  regard  had  to  rely 
on  calibration  measurements  on  a  homogeneous  me¬ 
dium  to  determine  these  critical  parameters.  Here 
we  show  that  we  can  use  a  simple  least-squares  min¬ 
imization  scheme  to  determine  these  parameters  ac¬ 
curately.  — 

To  illustrate  this  scheme,  we  show  how  one  can 
determine  the  BC  coefficient  a.  (The  source  term 
and  initial  estimates  of  the  optical  properties  can  be 
calculated  in  the  same  manner.)  We  compute  the  X2 
value  as  a  function  of  a. 


^  =  E  [$i(m>  -  <f>>>]2, 

i-1 

where  AT  is  the  number  of  boundary  measurements, 
<t>i m  is  the  measured  photon  density  from  a  given 
experimental  inhomogeneous  medium,  and  $/c)  is 
the  computed  photon  density  from  a  numerical  sim¬ 
ulation  of  a  homogeneous  medium  with  the  same 
geometry  as  the  experimental  medium.  To  calculate 
,  we  need  .to  give  only  a  range  of  the  a  value  that 
covers  the  correct  BC  coefficient.  With  the  given  a 
each  time,  we  compute  from  Eq.  (3).  The  ratio¬ 
nale  of  this  scheme  is  based  on  the  argument  that  the 
minimum  ofZ2  corresponds  to  the  correct  value  of  a 
associated  with  the  experimental  inhomogeneous 
medium.  In  the  same  manner  the  source  term  and 
initial  estimates  of  the  optical  properties  can  also  be 
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Fig.  1.  Schematic  representation  of  the  experimental  setup. 


computed.  Note  that  this  scheme  is  suitable  for  all 
data  types  including  dc,  frequency-,  and  time-domain 
data. 

B.  Experimental  Procedures  and  Parameters 

Our  experimental  setup,  shown  in  Fig.  1,  is  an  auto-~ 
mated  multichannel  frequency-domain  system  that 
has  been  described  in  detml  elsewhere.25  With  this 
system  a  radio-frequency  intensity-modulated  near- 
infrared  beam  at  785  run  is  sent  to  the  phantom  by  16 
fiber  optic  bundles  coupled  with  a  high-precision 
moving  stage.  The  diffused  light  is  received  by  an¬ 
other  16  channel  fiber  optic  bundles  and  delivered  to 
a  photomultiplier  tube  (PMT).  A  second  PMT  is 
used  to  record  the  reference  signal.  These  PMTs 
are  supplied  at  a  radio-frequency  modulated  current 
with  a  0 . 1-1-kHz  shift.  The  intermediary  frequency 
signal  obtained  from  the  PMT’s  is  processed  with  a 
National  Instruments  board.  To  increase  the  dy¬ 
namic  range  of  intensities  that  the  PMT  can  detect, 
an  automated  filter  wheel  with  precalibratedTheutral 
density  filters  is  added-fo  the  system.  For  every 
source  position,  16  measurements  for  each_detector 
were  made,  taking  alternately  100-ms  samples  for 
— sample  and  reference  signals.  The  dc,  ac,  and  phase 
shift  between  reference  and  sample  signals  were  ob¬ 
tained  by  fast  Fourier  transform  Labview  routines. 
The  total  data  collection  time  for  256  measurements 
was  8  min.  We  used  only  measured  dc  data  to  re¬ 
construct  the  absorption  and  scattering  images  in 
this  study. 

In  our  experiments  we  used  a  50-mm-diameter  cy¬ 
lindrical  solid  phantom  (1%  Intralipid  +  India  ink  + 
Agar)  as  the  background  medium.  In  the  cases-ex- 
amined  the  absorption  coefficient  pa  of  the  back¬ 
ground  varied  from  0.0005/mm  to  0.005/mm,  while 
the  reduced  scattering  coefficient  \l$'  was  always  1.0/ 
mm.  Both  single-  and  two-target  configurations 
were  used.  In  the  single-target  cases,  target  sizes  of 
15.0,  7.0,  and  4.0  mm  were  examined.  In  the  two- 
target  case,  one  target  6.0  mm  in  diameter  and  one 
8.0  mm  in  diameter  were  embedded  in. the  back¬ 
ground.  The  jxc  of  the  target  varied  from  0.015/mm 
to  0.050/mm;  the  pus'  of  the  target  was  1.0,  1.25,  or 
1.50/mm.  The  two-dimensional  finite  element 
mesh  used  had  249  nodes  and  448  elements  for  both 
forward  and  inverse  solutions.  Note  that  here  we  do 


Fig.  2.  (a)  Three-dimensional  plot  of  the  X2  error  versus  the  BC 
coefficient  a  and  absorption  coefficient  \ia.  (b)  X2  error  versus  the 
BC  coefficient  a  at  a  different  absorption  coefficient,  (c)  JX2  error 
versus  source  strength,  (d)  Reconstructed  p.a  image,  (e)  Target/ 
background  geometry  and  reconstructed  profile  of  p.a  along  the 
cut-line,  AB.  In  this  case  a  15-mm  single  target  with  .  - 
1.0/mm  and  [ia  =  0.020/mm  was  embedded  in  a  background  with 
=  1.0/mm  and  (±^-=  0.0005/mm.  _ 


not  need  to  use  the  dual-mesh  scheme  as  used- be¬ 
fore20-21  to  obtain  quantitative  reconstruction  by  dc 
data.  All  the  reconstructions  were  results  of  20  it¬ 
erations,  after  which  no  noticeable  improvement  was 
observed.  The  computations  were  performed  on  a 
400-MHz  Pentium  II  personal  computer. 

3.  Results  and  Discussion 

A  series  of  experiments  was  conducted  to  evaluate 
the  capabilities  of  our  improved  reconstruction  algo¬ 
rithm  to  detect  heterogeneities  in  the  presence  of 
single  and  multiple  targets  with  different  sizes  and 
contrasts  relative  to  the  background.  Both  color- 
scale  images  reconstructed  and  quantitative  plots  of 
one-dimensional  profiles  of  the  exact  and  recovered 
property  distributions  contained  in  these  images  are 
presented  in  this  section.  _ 

Through  a  simple  single-target  experiment,  in  Fig. 
2  we  present  detailed  results  from  data  preprocess¬ 
ing/optimization.  In  this  case  the  target  was  an  ab¬ 
sorber  with  only  jxa  =  0.02/mm.  Figure  2(a)  shows 
a  three-dimensional  plot  of  the  A2  error  versus  the  B 
coefficient  a  and  the  initial  estimate  value  of  p-o* 
(Because  we  used  the  same  p./  for  the  target  an 
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background  in  this  study,  we  easily  found  its  value  to 
be  0.99/mm  by  using  this  preprocessing  procedure.) 
As  can  be  seen,  a  minimum  value  of  the  X2  error  is 
clearly  reached  at  a  =  0.55  and  =  0.002/mm.  To 
examine  this  three-dimensional  plot  closely,  one¬ 
dimensional  profiles  of  the  X2  error  as  a  function  of  a 
are  provided  at  a  different  value  of  \ia  in  Fig.  2(b). 
Again,  the  minimum  is  confirmed.  Here  when  <J> 
(the  photon  density  for  a  numerical  simulation  homo¬ 
geneous  medium  with  the  same  geometiy  as  the  ex¬ 
perimental  medium)  was  computed,  the  source  term 
in  Eq.  (1)  was  already  normalized  with  the  highest 
photon  density  at  the  boundary.  .With  the  optimal 
values  of  a  and  |xa  determined,  the  optimal  source 
strength  SQ  can  now  be  obtained  by  computing 
with  an  unnormalized  source  term  in  Eq.  (1).  Figure 
2(c)  gives  the  result  from  this  calculation  in  which  the 
minimum  of  the  X2  error  is  reached  at  S0  =  49,500 
photons /mm3 /s.  Using  these  values  of  critical  pa¬ 
rameters  determined  from  the  data  preprocessing/ 
optimization  scheme,  we  reconstructed  the 
absorption  image  as  shown  in  Fig.  2(d).  The  absorp¬ 
tion  profile  along  the  cut  line  AB  is  displayed  in  Fig. 
3(e).  It  is  clear  that-the  recovered  image  is  quanti¬ 
tative  in  terms  of  the  optical  property  value,  target 
location,  and  size. 

Figure  3  shows  the  results  from  another  two  sets  of 
experiments  with  the  same  target  configuration  as 
the  case  above  but  with  different  values  of  \±a 
[0.015/mm  in  Figs.  3(a)  and  3(b)  and  0.035/mm  in_ 
Figs.  3(c)  and  3(d)].  As  can  be  seen  in  Figs.  3(a)  and 
3(c),  the  minimum  error  was  reached  at  a  =  0.553 
and  fxa  =  0.004/mm  and  at  a  =  0.540  and  = 
0.0075/mm,  respectively.  These  values  of  a  and  jxa 
were  then  used  to  obtain  the  optimaTsource  strength 
values.  The  reconstructed  optical-images  and  the 
target  absorption  profiles  are  shown  in  Figs.  3(b)  and 
3(d).  As  shown,  the  images  are  quantitatively  re¬ 
constructed  in  terms  of  the  localization,  size,  and  op¬ 
tical  properties  of  the  target. 

The  capability  of  our  improved  algorithm  to  recon¬ 
struct  a  small  target  with  both  \La  and  jjls'  different 
from  the  background  values  is  examined  in  Fig.  4. 
hi  this  experiment  a  4.0-mm-diameter  target  with  jx/ 

~  h25/mm  and  fia  =  0.050/mm  was  embedded  in 
the  background  medium  with  \i$'  =  1.0/mm  and  jxa  = 
0.0005/mm.  Both  reconstructed  images  and  prop- 
erty  profiles  for  both  fx^Land  |xa  are  presented  in  Fig. 

4-  Note  that  the  p,s'  image  can  be  recovered  quan¬ 
titatively  even  though  its  contrast  with  the  back¬ 
ground  is  low  (1.25:1).  The  boundary  artifacts 
appearing  in  the  scattering  image  could  be  partially 
ue  to  the  a  small-target  perturbation  that  in  this 
°ase  led  to  overall  low  signal  sensitivity  for  the  re¬ 
construction.  The  inherent  ill-conditioned  nature 
6  inverse  Problem  that  we  have  here  makes  it 
Ghiicult  to  recover  absorption  and  scattering  images 
^rnultaneously.  It  appears  to  be  generally  more  dif- 
cult  for  the  scattering  image  to  be  reconstructed 
^han  the  absorption  image.20'23 
In  Fig.  5  are  displayed  the  reconstructed  images 
the  property  profiles  in  which  two  differently 


Fig.  3.  (a)  X2  error  versus  the  BC  coefficient  a  at  a  different 
absorption  coefficient,  (b)  Reconstructed  p.a  image  and  profile 
along  the  cut-line^AB.  (c)  X1  error  versus  the  BC  coefficient  a  at 
a  different  absorption  coefficient,  (d)  Reconstructed  p,a  image  and 
profile  along  the  cut-line,  AB.  The  background_optical  properties 
for  these  two  cases  are  p/  -  1.0/mm  and  pa  =  0.0005/mm.  The 
optical  properties  of  the  15-mm  target  are  p/  =  1.0/mm  and  pa  = 
0.015/mm  for  (a)  and  (b)  and  p,'  =  1.0/mm  and  pa  =  0.035/mm  for 
(c)  and  (d),  respectively. 


sized  (6-  and  8-mm)  targets  were  embedded  in  the 
background  medium  (p.s#  =  1.0/mm  and  ^7  = 
0.0005/mm).  The  two  largets  have  the  same  ^  as 
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Fig.  4.  Reconstructed  p,'  and  pa  images  and  their  profiles  along 
the  cut-line,  AB,  for  a  4- mm  target  that  was  embedded  in  the 
background  with  p,'  =  1.0/mm  and  p0  =  0.0005/mm.  The  target 
optical  properties  were  p/  =  1.25/mm  and  pa  =  0.040/mm. 
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Fig.  5.  Reconstructed  and  fxa  images  and  their  profiles  along 
the  cut-line,  AB,  for  two  differently  sized  (6-  and  8-mm)  targets 
that  were  embedded  in  the  background  with  ji,'  =  1.0/mm  and 
=  0.0005/mm.  The  left  6-mm  target  was  p.a'  =  1.0/mm  and 
=  0.030/mm,  and  the  right  8-mm  target  was  pts'  *  1.0/mm  and 
—  0.050/mm. 

the  background  but  different  fxa  (0.030/mm  for  the 
left  6-mm  target  and  0.050/mm  for  the  right  8-mm 
target).  As  can  be  seen,  both  targets  can  be  clearly 
detected  quantitatively  in  terms  of  the  optical  prop¬ 
erty  value,  size,  and  location  of  the  targets. 

The  last  experiment  was  a  configuration  in  which  a 
10-mm  target  with  both  and  \is'  different  from  the 
background  was  embedded  in  a  background  with  a 
larger  jxa  (0.005/mm)  than  that  of  the  background 
used  in  the  previous  cases.  The  background  had  \xs ' 
=  1.0/mm.  The  target  had  \ls'  =  1.5/mm.and  fia  = 
0  09,0 /mm.  The  reconstructed  images  and  property 
profiles  for  both  jxa  and  [x/  are-presented  in  Fig.  6. 
Again,  quantitative  reconstructions  are  in  evidence. 
Note  from  Figs.  4  and  6  that,  although  the  low- 
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Fig.  6.  Reconstructed  p./  and  fia  images  and  their  profiles  along 
the  cut-line,  AB,  for  a  10-mm  target  that  was  embedded  in  the 
background  with  p.s'  =  1.0/mm  and  [ia  =  0. 0050/mm.  The  target 
optical  properties  were  =  1.5/mm  and  p,a  =  0.020/mm. 


contrast  fxs'  images  can  be  quantitatively  recovered, 
it  is  difficult  for  the  jis'  profiles  to  fully  reach  their 
exact  values.  In  contrast,  the  p,a  profiles  can  fully 
reach  their  exact  values.  We  encountered  the  same 
situation  in  image  reconstructions  using  frequency- 
domain  data.19-21  Although  further  study  is  war¬ 
ranted  in  this  regard,  we  suspect  that  there  is  a 
mathematical  competition  between  \xsr  and  pa  due  to 
their  own  unique  positions  in  the  diffusion  equation. 
Improvement  in  image  quality  may  be  realized  if  we 
consider  this  competition  in  the  reconstruction  algo¬ 
rithm. 


4„  Conclusions 

We  have  demonstrated  that  the  quantitative  good 
reconstruction  of  absorption  and  scattering  images 
can  be  achieved  with  dc  data  only.  A  series  of  ex¬ 
periments  with  both  single-target  and  multitarget 
configurations  has  been  conducted  to  evaluate  our 
improved  reconstruction  algorithm.  The  data  pre¬ 
processing/optimization  scheme  developed  has  been 
examined  in  detail.  This  scheme  allows  us  to  deter¬ 
mine  accurately  all  the  critical  initial  parameters- 
needed  for  reconstruction:  the  BC  coefficient  re¬ 
lated  to  the  tissue-detector-air  interface,  the  source 
term  associated  with  the  model  mismatch,  and  the 
initial  estimates  that  are  critical  for  a  nonlinear 
optimization-based  algorithm.  In  a  practical  situa¬ 
tion,  such  as  an  m  vivo  measurement,  the  tissue- 
detector— air  interface  can  be  quite  complicated;  using 
Hie  previous  calibration  procedure19-21  to  obtain  an 
accurate  BC  coefficient  is  difficult.  In  fact  such  a 
calibration  procedure  requires  a  homogeneous  me¬ 
dium  measurement,  whichJmimpossible  in  a  clinical 
situation.  Similarly,  when  one  wants  to  reconstruct 
a  two-dimensional  image  from  measurements-per-  ^ 
formed  in  an  actual  three-dimensionalmitting,  the 
model  mismatch-always  exists.  This  study  shows 
that  the  proposed  simple  data  preprocessing  proce¬ 
dure  can  minimize  the  difference  owing  to  this  model 
mismatch.  Finally  this  procedure  also  reduces  com¬ 
putation  cost  significantly  because  it  allows  us  to  use 
a  coarse  single  finite  element  mesh  for  both  forward 
and  inverse  procedures. 
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.ABSTRACT:  In  this  article,  we  present  a  novel  approach  that  combineslhe  best  aspects  of  both 
endogenous  opticaTimaging  and  exogenous  fluorescent  lifetime  and  yield  imaging  for  breast  cancer 
detection  and  diagnosis.  Using  this  approach,  spatial  distributions  of  optical  properties,  fluorescence 
lifetime,  and  yielcTin  tissue  can  be  reconstructed  from  measured  excitation  and  emission  data  at  the 
surface  of  the  breast  tissue  by  regularized  inverse  algorithms.  We~Show  images  from  simulated  data, 
as  well  as  images  from  experimental  data  using  tissue-like  phantom  materials  in  the  laboratory. 

KEY  WORDS:  optical  imaging,  fluorescence  imaging,  frequency-domain,  breast  cancer. 


I.  INTRODUCTION  —  - 

X-ray  mammography  is  the  clinical  tool  currently  used  for  breast  cancer 
detection.  However,  mammography  is  ionizing  radiation  and  has  an  unacceptable 
false-negative  rate  for  patients  with  radiodense  breast  tissues.  These  patients 
represent  the  general  population  of  premenopausal  women  as  well  as  those  with 
fibrocystic  tissue  disease.  Yet,  cancer  in  younger  women  tends  to  be  more  virulent 
and  grow  faster.  Consequently,  there  is  a  critical  need  to  seek  alternatives  that 
could  overcome  the  problems  associated  with  X-ray  mammography  and  that  could 
become  an  excellent  adjunct  or  competitive  tool  to  conventional  X-ray  mammog¬ 
raphy.  Nonionizing,  noninvasive  optical  imaging  has  great  potential  to  become  one 
of  the  most  promising  alternative  strategies  for  breast  cancer  detection.  This  is 
underscored  by  the  fact  that  in  the  last  few  years  significant  progress  has  been  made 
in  understanding  the  fundamental  nature  of  light  propagation  in  tissue,  which  has 
led  to  the  emergence  of  several  advanced  techniques  for  transmission  of  light 
through  tissue  over  distances  relevant  to  breast  imaging.1  Considerable  effort  has 
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been  made  toward  the  realization  of  a  new  class  of  breast  cancer  detection  methods 

using  these  emerging  optical  techniques. 

Both  direct  imaging  and  image  reconstructions  using  multiply  scattered  light 
have  been  demonstrated,  largely  from  simulations  and  tissue  phantom  experi¬ 
ments.  The  majority  of  these  optical  imaging  techniques  employ  methodologies 
that  are  based  on  differences  in  the  absorption  and/or  scattering  properties 
between  the  embedded  object  and  its  surroundings.  While  in  vitro  studies  show 
that  the  endogenous  optical  contrast  between  the  tumor  and  normal  tissues  can 
be  quite  low,2  in  vivo  experiments  indicate  that  this  endogenous  optical  contrast 
can  be  as  high  as  fourfold.3  Nonetheless,  development  of  methods  for  improving 
both  the  sensitivity  and  specificity  for  detection  is  necessary.  Fluorescence 
spectroscopy/imaging  studies  of  human  tissue  suggest  that  a  variety  oflesions 
show  distinct  fluorescence  spectra  compared  with  those  of  normal  tissue.  It  has 
been  shown  that  exogenous  dyes  may  offer  the  best  contrast  for  optical  imag¬ 
ing.8-9  Use  of  exogenous  agents  would  provide  fluorescent  markers  that  could 
serve  to  detect  embedded  tumors  in  the  breast.  In  particular,  the  ability  to  monitor 
the  fluorescent  yield  and  lifetime  may  also  provide  biochemical  specificity  if  the 
fluorophore  is  sensitive  to  a  specific  metabolite,_such  as-erxygen,  glucose,  etc.10 
This  has  intrinsic  merit  that  extends  beyond  the  structural  information  contained 
in  the  static  images  associated  with  X-ray  mammography,  thus  provides  another 
strong  rationale  for  its  complementary  role  in  the  detection/diagnosis  of  breast 
cancer.  Because  tumors  exhibit  distinctive  metabolite  alternations,  this  ability  to 
extract  tissue  functional  information  may  offer  a  unique  way  for  noninvasively 
differentiating  between  benign  and  malignant  tumors. 

In  this  paper,  we  present  a  novel  approach  that^combines  the  best  aspects  of 
both  endogenous  optical  imaging  and  exogenous  fluorescent  lifetime  and  yield 
imaging  for  breast  cancer  detection  and  diagnosis.  Using  this  approach,  spatial 
distributions  of  optical  properties,  fluorescence  lifetime,  and  yield  in  tissue  can  be 
_  reconstructed  from  measured  excitation  and  emission  data  at  the  surface  of  the 
breast  tissue  by  regularized  inverse  algorithms.  We  reporHmages  from  simulated 
data,  as  well  as  images  from  experimental  data  using  tissue-like  phantom  materials 
in  the  laboratory/ 


II.  THEORETICAL  AND  EXPERIMENTAL  METHODS 

In  frequency-domain,  it  is  known  that  propagation  of  both  excitation  and 
fluorescent  emission  light  in  tissues  or  multiply  scattering  media  can  be  described 
by  the  following  coupled  diffusion  equations:11 


V-.  [Dx(r)V<S>x(r,  ©)]  - 
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where  Ox,m  is  the  photon  density  for  excitation  (subscript  x)  or  fluorescent  light 
(subscript  m),  Dx,m  is  the  diffusion  coefficient,  pa  is  the  absorption  coefficient 
due  to  contributions  from  both  nonfluorescing  chromophores  and  fluorescent  dye, 
pa  is  the  absorption  coefficient  for  the  excitation  light  due  to  contribution  from 
fluorescent  dye,  co  is  the  modulation  frequency,  c  is  the  velocity  of  light  in  the 
medium,  and  r|  and_t  are-the  fluorescent  quantum  yield  and  lifetime,  respectively. 
S(r,  co)  is  the  excitation  source  term  in  (1)  which  for  a  point  source  can  be  written 
as  S  =  5o8(r  -  ro),  where  So  is  the  source  strength  and  8(r  -  ro)  is  the  Dirac-delta 
function  for  a  source  at  ro.  Note  that  a  single-exponential  fluorescence  decay  has 
been  assumed  in  the  source  term  for  fluorescent  light  (right-hand  side  of  (2)); 
multiexponential  time  decay  can  be  handled  by  a  simple  extension.  The  diffusion 
coefficient  can  be  written  as  Dx,m(r)  =  l/3(|i  (r)  +  u'  (r)],  where  u.'  (r)  is 

1  1  .  ax,m  ,  “r.m  — 

the  reduced  scattering  ^coefficient.  For  known  optical  properties  and-  fluorescent 
lifetime  and  yield,  Eqs.  (1)  and  (2)  become  standard  boundary  value  problems  for  the 
spatially  varying  photon  densities  of  excitation  and  emission  light  subject  to  appro- 
-  priate  boundary  conditions  (BCs).  We  use  the  non-zero  photon  density  BCs  or  type 
III  BCs  in  this  article,  which  are  proved  to  be  the  most  accurate  ones:12 
~DxmV<&xm  ■  n  =  aOx  m,  where  n  is  the  unit  normal  vector  for  the  boundary 
~  surface  and  a  is  a  coefficient  that  isrelated  to  the  internal  reflection  at  the  boundary. 

In  combined  optical  and  fluorescence  imaging,  the  goal  is  to  recover  all  distri¬ 
butions, ^including  Dx,-m,  |ia  ,  x,  and  r\.  Our  image  reconstruction  algorithms  are 
centered  on  the  finite  element  method.  The  image  formation  algorithms  cast  the 
inverse  problem  associated  with  determining  the  optical  and  fluorescent  characters 
of  the  medium  (i.e.,  breast  tissue)  as  a  nonlinear  parameter  estimation  where  known 
excitations  accompanied  with  measured  data  are  used  to  find  a  “best  fit”  of  the  tissue 
optical  and  fluorescent  parameters  needed  to  reproduce  the  known  information.  The 
details  of  our  reconstruction  algorithms  have  been  described  in  References  1 2  and  1 3 . 

The  experimental  system  used  for  the  combined  optical  and  fluorescence 
imaging  is  a  frequency-domain  optical  system  that  has  been  described  in  detail  in 
Reference  12.  Usingjhe  frequency-domain  imaging  system,  the  phase-shift  and 
intensity  for  excitation  and  fluorescent  light  can  be  measured. 


III.  RESULTS  AND  CONCLUSIONS 

For  simulations,  Figure  1  shows  the  lifetime  and  yield  images  reconstructed 
under  conditions  of  no  noise  and  with  5%  added  noise  for  the  “measured”  intensity 
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FIGURE  1.  Simulated  simultaneous  reconstruction  of  both  fluorescent  lifetime  (x,  ns)  and 
yield  (ri,  dimensionless)  under  different  noise  condition,  (a)  x  reconstruction  with  no  noise 
added;  (b)  x  reconstruction  with  5%  random  noise  added;  (c)  n  reconstruction  with  no  noise 
added;  (d)  n  reconstruction  with  5%  random  noise  added.  In  these  simulations,  a  15-mm- 
diameter  object  is  embedded  in  a  43-mm-diameter  background  region. 


and  phase  shift  of  both  excitation  and  emission  light.  As  can  be  seen,  the  images 
formed  are  qualitatively,  correct,  even  for  those  with  a  5%  noise  level.  For  experi¬ 
ments,  Figure  2  presents  optical  property  images  reconstructed  using  tissue-like 
phantom  materials  that  simulate  soft  tissues.  Here  we  show  that  a  4-mm-size  object 
can  be  clearly  recovered  using  our  reconstruction  approach. 

In  summary,  wertiave  demonstrated  optical  and  fluorescent  image  reconstruc¬ 
tions  using  simulated  and  experimental  data.  Our  results  show  that  it  is  possible  to 
detect  breast  tumors  using  our  approach  described  here.  At  Clemson  we  are  currently 
developing  a  clinical  prototype  imaging  system  for  human  subject  studies.  _ 


ACKNOWLEDGMENTS  — 

This  work  was  supported  in  part  by  the  National  Institutes  of  Health  (NIH) 
(R55  CA78334),  Oak  Ridge  Associated  Universities,  and  the  Greenville  Hospital 
System/Clemson  University  Biomedical  Cooperative. 


FIGURE  2.  Optical  property  images  reconstructed  from  experimental  data  using  tissue¬ 
like  phantom  materials  (Intralipid  mixed  with  India  ink).  The  diameter  of  the  background 
region  is  86  mm.  The  data  used  were  collected  from  the  frequency-domain  imaging  system 
described  in  Reference  12. 
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ABSTRACT"  _ 

The  reconstruction  of  fluorescence  lifetime  distributions  in  heterogeneous  turbid_media  and  tumor-bearing  animals  are 
experimentally  demonstrated  by  frequency-domain  measurements.  A  set  of  coupled  diffusion  equations  are  used  to  describe 
the  propagation  of  excitation  and  fluorescent  emission  light  in  multiply  scattering  media.  A  finite  element  based  reconstruction 
algorithm  combined  with  Marquardt  and  Tikhonov  regularization  methods  are  used  to  obtain  the  fluorescence  images.  The 
experimental  set-up  is"  an  automatic  multi-channel  frequency-domain  system.  16  sources  and-  16  detectors  are  used. 
Experiments  are  performed  using  indocyanine  green  (ICG)  and  3,3’-diethylthiatricarbocyanine  iodide  (DTTCI)  in  tissue-like 
phantoms  of  both  single-  and  multi-target  configurations  with  considerations  of  perfect  and  imperfect  uptake  of  fluorescence 
dyes  in  the  scattering  media.  ICG  are  used  in  tumor-bearing  animal  studies.  Our  results  show  that  the  fluorescence  lifetime 
image  of  the  heterogeneities  within  a  circular  surrounding  medium  and  in-vivo  tissue  can  be  reconstructed  successfully. 

1.  INTRODUCTION 

Recently  the  idea  of  reconstruction  life-based  fluorescence  diffusion  tomography  has  been  proposed  and  developed  for  clinical 
applications  such  as  breast  cancer  detection  and  tissue  functional  mapping.1'6  This  new  imaging  approach  relies  on  the  fact 
that  the  lifetime  of  fluorophores  in  tissue  can  potentially  provide  tissue  functional  information  such  as  tissue  oxygenation,  pH, 
and  glucose.7  It  is  also  based  on  the  fact  that  the  fluorophore  will  preferentially  accumulate  in  tumors,  hence  providing  the  best 
sensitivity  for  cancer  detection.8  In  this  type  of  indirect  imaging  method,  an  effective  reconstruction  algorithm  is  crucial  which 
allows  for  the  formation  of  a  spatial  map  of  the  lifetime  of  heterogeneous  fluorophore  distributions.  To  date  fluorescence 
image  reconstructions  are  largely  limited  to  simulated  data,  based  on  linear  and  nonlinear  algorithms.1'6  Chang  et  al.2 
attempted  to  obtain  a  qualitatively  good  image  of  the  fluorophore  concentration  only  from  dc  measurements.  In  this  paper  we 


present  successful  reconstruction  of  lifetime  images  from  frequency-domain  measurements  using  indocyanine  green  (ICG) 
dye  and  3,3’-diethylthiatricarbocyanine  iodide  (DTTCI)  in  tissue-like  phantoms  of  single-  and  multi-target  configurations  with 
considerations  of  perfect  and  imperfect  uptake  of  fluorescence  dyes  in  the  scattering  media.  And  ICG  are  used  in  tumor¬ 
bearing  animal  studies. 


2.  RECONSTRUCTION  ALGORITHM 

In  frequency-domain,  it  is  known  that  propagation  of  both  excitation  and  fluorescent  emission  light  in  tissues  or  turbid  media 
can  be  described  by  the  following  coupled  diffusion  equations:2'6 
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where  Oxm  is  the  photon  density- for  excitation  (subscript  x)  or  fluorescent  fight  (subscript  m),  Dx  m  is  the  diffusion 
coefficient,  (Ia  is  the  absorption  coefficient  due  to  contributions  from  both  non-fluorescing  chromophores  and  fluorescent 
dye,  (Ia  is  the  absorption  coefficient  for  the  excitation  light  due  to  contribution  fromTTuorescent  dye,  CO  is  the  modulation 
frequency,  c  is  the  velocity  of  light  in  the  medium,  and  T|  and  X  are  the  fluorescent  quantum  yield  and  lifetime,  respectively. 
S(r,CO)  is  the  excitation  source  term  in  (1).  Note  that  a  single-exponential  fluorescence  decay  ha^ been  assumed  in  the  source 
term  for  fluorescent  light  (right-hand-side  of  (2));  multi-exponential  time  decay  can  be  handled  by  a  simple  extension.  The 
diffusion  coefficient  can  be  written  as  Dy  m  =  ^  ^  )  where  |LL^  ^  is  the  reduced  scattering  coefficient.  In  this 

study  a  point  excitation  source,  S  =  S05(r  -  ro  ) ,  and  Type  m  boundary  conditions  (BCs):  “ D VOx  m  *  n  =  0C<I>x  m ,  are 
used,  where  S0  is  the  source  strength,  n  is  the  unit  normal  vector  for  the  boundary  surface  and  CX  is  a  coefficient  that  is 
related  to  the  internal  reflection  at  the  boundary.  Both  S0  and  a  are  determined  directly  from  the  heterogeneous  measurements 
by  use  of  a  data  pre-processing  scheme  that  is  reported  elsewhere.9 


In  fluorescence  lifetime  tomography,  the  goal  is  to  recover  all  distributions  including  Dxjn,  T  and  p.  We  are  focusing  on 
the  reconstruction  of  t.  Other  parameters  can  be  easily  dealt  with  in  the  same  way.  Making  use  of  finite  element  discretization, 
we  can  obtain  two  matrix  equations  for  Eqs.  (l)-(2)  and  realize  other  derived  matrix  relations  through  differentiation: 


[AjOx}={bx) 

[AmK<i>mMbm} 


(3) 

(4) 


where  the  elements  of  matrices  [AJ  and  [Am]  are  respectively  (ax  m  )y  =<  ~D x  m  V\y j  *  V\j/j  —  (|Ias 


-)VjVi  >>and 


the  entries  in  column  vectors  {bx>m}  and  {Oxm}  are  (bx)j  =  -  <  Sxj/j  >  +0cX(<I)x  )j  j ¥jVi^s  » 

j=i 


K  N  L  L  M 

(bm)i  =-<lTik\|/k^ia  mI(Ox)jvj\(/i(l-icoXTI\|/1)/(l  +  co2(Zx1\K1)2)>+aX(^ >m)jfVjVids,  and 

k=!  j=l  1=1  1=1  j=l 

<£x>m  =  { (Ox>m ) , ,  (O  x  m  )2  -•  •  •  ( <5x,m  ) N }T ,  where  (  )  indicates  integration  over  the  problem  domain,  and  <&xm,  X  has 
been  expanded  as  the  sum  of  coefficients  multiplied  by  a  set  of  locally  spatially- varying  Lagrangian  basis  functions  \|/j,  \j/j 

and  \|/k.  j  expresses  integration  over  the  boundary  surface  where  Type  III  BCs  have  been  applied.  is  the  photon 

density  at  node  i,  N  is  the  node  number  of  a  finite  element  mesh  and  M  is  the  boundary  node  number. 

To  form  images  from  presumably  uniform  initial  guess,  we  us^the  following  equation  to  update  x  from  their  startingrvalues.  — 

3At  =  OLo)-0^  _(5) 

where  O  J ° }  and  are  the  observed  and  computed  fluorescent  data,  respectively.  3  is  the  Jacobian  matrix  consisting  of 
the  derivatives  of  ®m  with  respective  to  t  at  each  boundary  observation  node;  At  is  the  vector  that  expresses  perturbations 


But  as  the  approximation  of  principle,  the  influence  of  coherent  noise  and  system  noise,  the  "ill-condition'’  of  the  equation  (5) 
is  almost  inevitable.  We  use  Marquardt  and  Tikhonov  regularization  methods  to  stabilize  the  decomposition  of  the  square 
system  of  equations. 


(3r3  +  A/}lT  =  3r(<0)  -  Oif) 


where  /  is  the  identity  matrix  and  X  maybe  a  scalar  or  a  diagonal  matrix. 

When  a  weighted  sum  of  the  squared  difference  between  computed  and  measured  data  is  minimized,  a  lifetime  image  can  be 
obtained. 


3.  EXPERIMENTS  AND  METHODS 


« 


3.1  Experimental  Set-up 


The  basic  components  of  our  frequency-domain  experimental  set-up  is  shown  in  Fig.  1.  The  system  used  a  near-infrared  laser 
diode  to  provide  excitation.  The  laser  was  modulated  by  a  radio- frequency  sine  wave  of  frequency/  that  was  amplified  by  RF 
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Fig.  1  Experimental  Set-up 


amplifier.  The  modulated  light  was  collected  by  the  collimating  system  and  then  focused  down  into  the  2.47  mm  core  diameter 
fiber.  16  source  fibers  and  16  detect  fibers  are  combined  together  to  form  a  measuring  a  ring.  The  laser  beam  was  sent  to  the 
phantom  by  16  fiber  optic  bundles  coupled  with  a  high  precision  moving  stage.  The  diffused  radiation  was  received  by  another 
16  channel  fiber  optic  bundles  and  delivered  to  a  PMT.  A  second  PMT  was  used  to  record  the  reference  signal.  These  PMTs 
were  driven  by  another  amplified  synthesizer  signal  at  a  frequency  of  /  +0.0001  -0.001MHz  .  A  heterodyning  effect  was 
achieved,  which  produced  a  100Hz~1000Hz  signal.  Both  reference  signal  and  diffusion  signals  were  digitized  by  a  data 
acquisition  board  which  is  inserted  in  the  expanded  slot  of  the  PC  computer.  In  order  to  increase  the  system  dynamic  range,  an 
automated  filter  wheel  with  pre-calibrated  neutral  density  filters  was  added  to  the  system.  For  every  source  position,  diffusion 
and  reference  signals  were  taking  alternatively  for  every  detect  fiber.  Fluorescence  signals  were  obtained  through  an 
interference  filter,  which  was  placed  in  front  of  the  detection  PMT.  AC,  DC  intensities  and  phase  shift  between  reference  and 
diffusion  signals  were  obtained  using  FFT  Labview  routines.  For  each  measurement  position,  it  would  need  about  100  ms  to 


acquire  the  data.  For  total  256  measurements,  it  would  take  about  8  minutes  to  complete  the  data  collection.  In  animal 
experiments,  only  11  source  fibers  and  11  detector  fibers  were  used.  The  total  time  for  121  measurements  were  about  5 
minutes. 

3.2  Experimental  methods 

Solid  phantom10  was  used  to  mimic  the  human  tissue.  It  was  made  of  agar,  Intralipid,  black  ink  and  fluorescent  dyes.  The 
absorption  and  the  reduced  scattering  coefficients  are  linear  with  the  ink  and  Intralipid  concentrations,  respectively. 
Micromolar  ICG  and  DTTCI  dyes  were  added  in  the  tissue  phantom  to  provide  the  fluorescence  contrast.  Agar  is  used  to  make 
the  phantom  solid.  The  solid  phantom  is  consisted  of  a  cylindrical  background  and  a  cylindrical  heterogeneity. 

The  absorption  peaks  ofiCG  and  DTTCI  are  764  nm  and  780  nm,  respectively.  And  the  fluorescent  emissionj>eaks  of  them 

t 

are  803  nm  and  830  nm  respectively.  Their  lifetimes  in  water  are  0.56  ns  and  1.18  ns,  respectively.8  A  laser  diode  of 
wavelength  785  nm  was  used  to  excite  both  dyes,  and  the  emission  light  of  wavelength  at  830  nm  were  detected  for  both  of 
them  through  an  interference  filter  of  center  wavelength  830  nm  with  lOnm  bandwidth. _  ~ 

Wistar  furth  female  rats  with  MT/W9a-B  (mammary  carcinosarcoma  )  cell  line  developedln  their  right  udder  were  used  as  the 
animal  models.  Animals  were  studied  when  the  tumor  sizes  were  around  10-30  mnr'diameters.  Animals  were  initially 
anesthetizecLwith  2%  isoflurane.  The  femoral  vein  was  cannulated  for  the  administration  of  intravenous  drugs,  including  the 
dye  ICG.  Anesthesia  was  maintained  during  the  experiment.  The  animals  were  placed  in  a  holder  with  their  tumor  positions 
placed  intojhe  measuring  ring.  1.5  mg/Kg  bodyweight  ofiCG  dye  was  injected  into  the  animal.  Imaging  studies  were  carried 
out  after  2  minutes  of  the  dye  injection.  _  —  _ 

3.3  Experimental  Results 

3.3.1  Phantom  measurement  results 

Phantom  studies  were  carried  out  with  single  and  multi-targets  configurations.  Different  dye  concentration  contrast,  different 
absorption  coefficient  and  size  contrasts  between  the  background  and  the  heterogeneities  were  studied.  Perfect  (no  dye  was 
present  in  the  background)  and  imperfect  uptake  (background  with  dye)  of  fluorescence  dyes  in  the  scattering  media  also  were 
studied. 

Figs.  1-5  show  the  reconstructed  images  of  T  for  phantoms  with  one  or  two  off-centered  targets  .  Figs.  1-2  are  the  results  of 
perfect  uptake  of  the  fluorescent  dyes.  Figs.  3-5  are  the  results  of  imperfect  uptake  of  the  fluorescent  dyes,  with  lifetime  of  the 
target  is  shorter  or  longer  than  that  of  the  background.  As  can  be  seen,  all  the  targets  can  be  clearly  imaged.  The  targets' 
locations  and  sizes  are  also  correctly  resolved  (  exact  data  are  not  shown  ).  It  is  also  interesting  to  note  that  a  finite  value  of  X 
(larger  than  the  X  value  in  the  target)  for  the  background  region  is  always  reconstructed  in  the  perfect  uptake  cases,  although 
there  is  no  fluorophore  distribution  in  the  background. 


3.3.2  In  vivo  tissue  measurement  result 


Fig.  6  shows  a  reconstructed  image  of  T  for  a  tumor-  bearing  rat.  The  rat  has  a  body  size  of  33  mm  diameter  and  tumor  size 
of  10  mm  diameter.  The  weight  of  the  rat  is  155g.  The  reconstructed  lifetime  image  shows  that  the  lifetime  in  the  tumor  area  is 
longer  than  the  surrounding  normal  tissue.  ICG  dye  is  used  in  the  animal  studies.  The  reconstructed  lifetime  image  shows  the 
tumor  size  and  location  correctly  (exact  data  are  not  shown). 


Fig.  1  Reconstructed  lifetime  image  of  phantom.  The 
background  hadn’t  dye  inside  and  the  target  had  ICG  dye  of 
concentration  0.5  jiM.  The  diameter  of  the  background  and 
the  target  were  50.0  mm  and  12.8  mm,  respectively.  The 
reduced  scattering  coefficients  of  the  background  and  target 
were  the  same  as  0.5  mm'1.  The  absorption  of  the 
background  and  target  were  0.001  mm'1  and  0.005  mm'1  , 
respectively.  The  lifetime  of  ICG  in  water  is  0.56  ns. 


Fig.  2  Reconstructed  lifetime  image  of  phantom.  The 
background  hadn’t  dye  inside  and  the  targets  had  ICG 
dye  of  concentration  1.0  pM.  The  diameter  of  the 
background  and  the  target  were  50.0  mm  and  10.0 
mm,  respectively.  The  reduced  scattering  coefficients 
of  the  background  and  target  were  the  same  as 
0.5mm'1.  The  absorption  of  the  background  and 
target  were  0.001  mm'1  and  0.005  mm'1  , 
respectively.  The  lifetime  of  ICG  in  water  is  0.56  ns. 


Fig._  3  Reconstructed  -lifetime  image  of  phantom.  The 
background  had  ICG  dye  of  concentration  0.1  pM,  and  the 
target  had  DTTCI  dye  of  concentration  5.0  pM.  The  diameter 
of^he- background  and  the  target  were  50.0  mm  and  12.8  mmr 
respectively.  The  reduced  scattering  coefficients  of  the 
background  and  target  were  the  same  as  0.5  mm'1.  The 
absorption  of  the  background  and  target  were  0.001  mm’1  and 
0.005  mm’1 ,  respectively.  Theiifetimes  of  ICG  and  DTTCI  in 
water  are  0.56  ns  and  1.18  ns,  respectively. 


Fig.  5  Reconstructed  lifetime  image  of  phantom.  The 
background  had  ICG  dye  of  concentration  0.04  |iM,  and 
the  target  had  DTTCI  dye  of  concentration  1.21  pM.  The 
diameter  of  the  background  and  the  target  were  57.5  mm 
and  13.5  mm,  respectively.  The  reduced  scattering 
coefficients  of  the  background  and  target  were  the  same  as 
0.5  mm*1.  The  absorption  of  the  background  and  target 
were  0.005  mm’1  and  0.01  mm'1  ,  respectively.  The 
lifetimes  of  ICG  and  DTTCI  in  water  are  0.56  ns  and  1.18 
ns,  respectively. 


Fig.  4HReconstructed  lifetime  image  of  phantom.  The 
background  had  DTTCI  dye  of  concentration  0.1  pM,  and 
the  target  had  ICG  dye  of  concentration  3.0- pM.  The 
diameter  of  the  background  and  the  target- were  35.0-ram 
and  8.5  mm,  respectively.  The  reduced  scattering 
coefficients  of  the  background  and  target  were  the  same  as 
0.5  mm'1.  The  absorption  of  the  background  and  target 
were  0.001  mm'1  and  0.005  mm'.1  ,  respectively.  The 
lifetimes  of  ICG  and  DTTCI  in  water  are  0.56  ns  and  1.18 
ns,  respectively. 


Fig.  6  Reconstructed  lifetime  image  of  tumor-bearing 
rat.  The  diameters  of  the  tumor  and  the  rat's  body 
were  10mm  and  33mm,  respectively.  The  weight  of 
the  rat  was  155g.  The  lifetime  of  ICGjn  water  is 
0.56  ns. 


4  CONCLUSIONS 


Fluorescence  lifetime  images  reconstruction  in  turbid  media  and  in  vivo  tissues  were  studied  in  this  paper.  Both  single-  and 
multi-target  configurations  with  considerations  of  perfect  and  imperfect  uptake  of  fluorescence  dyes  in  scattering  media  were 
studied.  In-vivo  tissue  lifetime  image  reconstruction  were  studied  on  tumor-bearing  animals.  The  results  showed  that  the 
fluorescence  lifetime  of  the  heterogeneities  within  a  circular  surrounding  medium  and  in-vivo  tissue  can  be  reconstructed 
successfully  using  the  algorithm  and  frequency-domain  measurement  presented  in  this  paper. 
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